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Abstract 

Multipolar order in complex fluids is described by statistical correlations. This paper presents 
a novel dynamical approach, which accounts for microscopic effects on the order parameter space. 
y—{ Indeed, the order parameter field is replaced by a statistical distribution function that is carried by 

t— I the fluid flow. Inspired by Doi's model of colloidal suspensions, the present theory is derived from a 

hybrid moment closure for Yang-Mills Vlasov plasmas. This hybrid formulation is constructed under 
the assumption that inertial effects dominate over dissipative phenomena, so that the total energy is 
conserved. After presenting the basic geometric properties of the theory, the effect of Yang-Mills fields 
is considered and a direct application is presented to magnetized fluids with quadrupolar order (spin 
nematic phases). Hybrid models are also formulated for complex fluids with symmetry breaking. For 
\Q the special case of liquid crystals, the moment method can be applied to the hybrid formulation to 

study to the dynamics of cubatic phases. 
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1 Introduction 



1.1 Micromotion and its symmetry properties 

Complex fluids are fluid systems whose particles carry an internal structure. Interesting examples are 
found in the dynamics of ferromagnetic fluids and spin glasses [14, 23]. In this case, the particles 
carry a spin magnetic moment undergoing its own rotational motion (micromotion), usually under the 
influence of external magnetic fields. This micromotion exhibits a full rotational symmetry so that 
the corresponding dynamics is entirely expressed in terms of the spin angular momentum variable, a 
situation also occurring for the rotational motion of a free rigid body on the rotation Lie group SO(3). 
However, this fully symmetric configuration differs from many situations occurring in condensed matter 
physics, which exhibit the symmetry breaking phenomenon. The most celebrated example is provided 
by the dynamics of nematic molecules, which are endowed with a special orientation identified with an 
unsigned unit vector n ~ — n (director). The presence of this special direction is responsible for symmetry 
breaking: the micromotion is no more symmetric under all possible SO(3) rotations in physical space, 
but rather the molecule exhibits only an 0(2) residual symmetry. Then, the dynamics of the nematic 
molecule is described in terms of two variables: its director and the associated angular momentum. 

1.2 Correlation effects: complex media with multipolar order 

In common continuum descriptions of complex fluids, the effects of micromotion are introduced through 
a macroscopic continuum field taking values in the space M of internal dynamical variables (e.g. the spin 
or the director). This space is usually known as order parameter space. More particularly, one introduces 
a map <f> : R 3 — > M (order parameter field) from physical space M 3 to the order parameter space M. This 
quantity is assumed to arise from an averaging process involving an appropriate statistical distribution 
function p G C°°(R 3 x M). While this averaging process has been widely used in equilibrium statistical 
approaches since Onsager's work [52], its use in the study of dynamical aspects has always been difficult. 
As a result, microscopic dynamical properties of orientational order are neglected in many complex fluid 
models. 

On the other hand, an increasing number of complex media are being discovered in which statistical 
correlation effects determine the particle interaction. For example, magnets exhibiting quadrupolar order 
(spin correlations) known as "spin nematics" have been investigated since the early 80's [ ] and their 
dynamical description remains an open question in soft matter modeling [ ]. These systems possess 
an experimental realization, which was reported in [43]. Recently, investigations have begun that study 
magnetic media with multipolar order higher than two: in [ ] a phase diagram was presented to identify 
quadrupolar, octupolar and hexadecupolar ferromagnetic phases. Higher-order correlation effects also 
play some role in the dynamics of liquid crystals. Indeed, when molecules possess a special shape (e.g. 
hard-cut spheres), the emergence of a "cubatic phase" in liquid crystals [58, 11] has been observed in 
computer simulations for about 20 years. This phase owes its name to the cubic orientational order 
(third-order director correlations) of its molecular interaction. The word 'cubatic' is used to distinguish 
this phase from cubic phases, which possess also translational order. Cubatic phases have been shown to 
emerge for a variety of molecular shapes, see e.g. [2] for how cubatic phases emerge in cubelike molecules. 

The effects of multipolar order are usually studied by introducing a tensor order parameter field, 
which encodes the correlation effects. For example, this is the approach adopted for cubatic liquid 
crystals [11] and magnets with quadrupolar order [38]. The question whether such a description is 
satisfactory may require complicated kinetic theory arguments based on BBGKY-like hierarchies. This 
is precisely the case for Doi's celebrated description of colloidal solutions [10]. Doi's theory aims to 
couple the dynamics of an ordinary fluid with a Smoluchowski kinetic equation for the dynamics of 
colloidal particles with orientation. In this way, the fluid correlation effects are encoded as usual in the 
equation of state, while the statistical information on the colloidal particles is carried along by the fluid 
through the Smoluchowski distribution function. This description stands as the main inspiration for 
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the present work, whose purpose is to introduce a unified systematic framework for the hybrid fluid- 
kinetic description of perfect complex fluids, by invoking well established Lie symmetry principles. This 
framework can be used to model the fluid stress tensor, which remains an outstanding open question 
for these systems. An important contribution is due to Constantin [8], who derived a general form of 
the stress tensor from energy principle considerations. However, contrarily to Doi theory, the present 
treatment considers inertial effects that dominate over dissipation, so that the total energy is conserved. 

Besides Doi's approach [10], other attempts to account for correlation effects include mesoscopic 
theories of liquid crystals [3, 4]. These theories involve a set of fluid balance equations on the Cartesian 
product K 3 x M of physical space and order parameter space. However, these equations are enormously 
difficult to handle since each different fluid variable (e.g. density and momentum) is defined on a space 
with a high number of dimensions. 

Additional considerations are needed when orientational defects (i.e. disclinations) are present in the 
system. In this case, the interaction between fluid, order parameter variables and defect dynamics is of 
central interest in modern soft matter physics. To the author's knowledge, no microscopic description 
of this interaction is available in the literature. This description is one of the aims of the present paper. 

Before introducing the hybrid approach to multipolar interactions, the next section reviews the Yang- 
Mills- Vlasov system [21] governing, e.g. quark-gluon plasmas. Its fluid closure (in particular the MHD 
limit [ ]) provides the geometric basis of many complex fluid models that inspired the present work. 



1.3 Yang-Mills Vlasov plasmas: a review 



The kinetic theory of Yang-Mills plasmas finds its roots in the definition of a distribution function 
/(x, p, a, t), where (x, p) G T*IR 3 = M 6 are the usual position-momentum coordinates, while a G o* 
is the single-particle Yang-Mills charge belonging to the dual space of some Lie algebra o. When the 
standard Klimontovich method of kinetic theory [41] is applied in this setting, one finds the mean- field 
kinetic equation (Yang-Mills Vlasov, or YM Vlasov) [21, 50] 
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where {•, •} denotes the canonical Poisson bracket on 



and [•, •] denotes the Lie bracket on o, while 



(•, •) is the pairing between o and its dual. This equation is coupled to the Yang- Mills field equations to 
form the Yang-Mills Vlasov system, with Hamiltonian 
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where (A,E) G ^(M^o) x f2 2 (M 3 ,o*) are the YM field variables consisting of the vector potential A 
and the electric field part E. Here J7 1 (M 3 ,o) and f2 2 (M 3 ,o*) are o-valued one-forms and o*-valued two- 
forms respectively (f2 n (M) denotes a generic n-form on the manifold M), while d A ■ = V • — [A, •] is 
the covariant differential (sometimes also denoted by V A ). Also, Aq G C°°(IR 3 , o) is the o-valued scalar 



potential so that setting 5H/5Aq 
The YM field equations are 



produces the nonholonomic constraint div E+ J afd 3 pda = 0. 



d t E- ad* Ao E + J af (p - (a, A)) da = curl A B , d t A 



-E- 



-d A A 



(3) 



where B = d A A is the magnetic component of the YM field and the suffix A on differential operators 
always denotes covariant differentiation [21]. In particular, one has div A E = divE + ad^. E{ and 
curl A B = curlB — [^4j,i?fc], where summation over repeated indexes is intended. The above YM 



3 



field equations are coupled to the Vlasov equation 
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where we recognize that the Klimontovich particle solution [41, 21] 

/(x, p, a, t) = 5(x - Q(t)) <5(p - P(t)) <5(a - T(t)) 

consistently returns Wong's equations for a single charge moving in a Yang-Mills field [49, 50]. 

The equations of chromohydro dynamics [21] are obtained from the equation hierarchy for the mo- 
ments 

X n>k (x,t) = J p n a k f{x,p,a,t) dpda (5) 



by making use of the cold-plasma closure 
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Here, the zero-th moment Xq^ is evidently the particle density, while Xi o and Xq^i are the fluid mo- 
mentum and the Yang-Mills charge density respectively. The moments X nk possess a Poisson bracket 
[22], which is obtained by using the chain rule 
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n,k n,k 

(here _l denotes tensor contraction) in the following Poisson structure of the YM Vlasov equation [2] 



{F,K}(f)=j /|^,^|d 3 xdpd C T + J f(a 



d 5F d 5K 



d x dp da . 



(8) 



Then, retaining only the terms involving D = Xq } o, G = Xq \ and m = X^o produces the following 
equations of chromohydrodynamics [21] 
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In this paper, Lie derivatives along the vector field SH/5m are either denoted by £sH/Sm or expressed 
explicitly, depending on convenience. The YM Vlasov Hamiltonian transforms into 

H =\ j l m - (G, A)| 2 d 3 x + J D U(D) d 3 x 

-J {G,A )d 3 x+^ J (\d A A\ 2 + \E\ 2 ^jd 3 x- J (E,d A A )d 3 x (10) 

where the internal energy D U(D) is inserted to account for pressure effects. This yields the following 
equation for the velocity u = D^ 1 (m — (G, A)): 



du 
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where p = D 2 



dU/dD. This equation is coupled to the Yang-Mills field equations 
d t E- adV E + Gu = curl A B , d t A = -E-d A A 



(12) 



and the transport equations for the fluid density D and the charge density G: 



d t D + div(L> u) = , 



d t G + adV G + div A 



(Gu) = 



(13) 



Here the "hydrodynamic gauge" Aq + u • A = is adopted. 

An important result concerning the CHD equations is that, in their MHD limit [3 L], their Hamiltonian 
structure is identical to the equations for spin glasses, with the only difference residing in the explicit 
expression of the Hamiltonian [30]. Without entering the treatment of Yang-Mills MHD, the next 
section introduces the hybrid closure of the Yang-Mills- Vlasov distribution function. This closure is the 
basis for all the complex fluid models presented later. 

1.4 Hybrid models 

This paper provides a general setting to include microscopic effects in the micromotion of complex fluids 
by the introduction of a new moment closure of the Yang-Mills Vlasov kinetic equation. This closure 
is similar to the cold plasma closure (6) that arises in the dynamics of quark-gluon plasmas in neutron 
stars. The motivation for this approach is the well known correspondence between many soft matter 
models and the fluid theory underlying Yang-Mills collisionless plasmas, which was reviewed in the 
previous section. In particular the equations (9) of chromohydrodynamics have been shown to apply 
also to complex fluids such as spin glasses and liquid crystals [24, 15]. 

This paper develops new theories of perfect complex fluids in which the transported order parameter 
field </> : M 3 — >• M is replaced by a position- dependent probability distribution function p £ C°°(M 3 x M) 
on the order parameter space M. This statistical quantity arises very naturally in the context of kinetic 
equations of the type Vlasov- Yang-Mills, for which M is a Lie group often called order parameter group 
O [24, 17]. When particle dynamics is O-invariant, simple Lie-Poisson reduction (i.e. T*OjO = o*) 
takes the dynamics on the cotangent bundle T*0 to the reduced dynamics on the dual Lie algebra o*. 
Then, a new fluid model arises from a moment closure of the Vlasov distribution function, which is 
defined on T*(R 3 x 0)/0 ~ M 6 x o* (see previous section). In particular, if / £ C°°(IR 6 x o*) is the 
Vlasov distribution, one can introduce the cold-plasma-like closure 



One notices that the distribution p(x, a,t) = ff dp belongs to a moment hierarchy that is different 
from the one containing the other two moments m(x, t) = Jp/dpdcr and D(x,t) = J f dp da, since 
p does not involve integration over a 6 o*. From a mathematical point of view, the variable p is a 
probability distribution function on o* that is attached to each fluid particle moving in M 3 . Hence, we 
write p G Den(IR 3 x o*), where Den denotes a space of probability densities. Consequently, one can 
consider p as a map belonging to the dual space of maps IR 3 — > C°°(o*). On the other hand, C°°(o*) is 
identified with the Poisson algebra of O-invariant Hamiltonian functions on T*0 (here also denoted by 
C°°(T*0)o). Thus, we shall see that from a mathematical point of view, the introduction of the variable 
p amounts to introducing the infinite-dimensional order parameter space C°°(o*), so that the geometric 
properties of complex fluid dynamics [24, 15] can be naturally extended to the case of infinite-dimensional 
order parameters. 

The special feature of this approach is that the resulting fluid model naturally accounts for the 
statistical information on the micromotion. This is a statistical distribution which is dragged along the 
fluid particle trajectories, as it also happens for ordinary density variables in standard hydrodynamic 
models. Because of the coexistence of both the statistical distribution on M and the ordinary fluid 
variables such as momentum and mass density, we call these models hybrid models. 
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Figure 1: Illustrative comparison of the conventional and the hybrid approaches with order parameter 
space M = S 2 . While the conventional approach associates a director n(x) to each fluid parcel at position 
x, the hybrid approach associates the whole probability density p(x, n). Probability densities of the same 
type also appear in Doi's theory of polymeric fluids [10] with negligible inertial effects. 

Plan of the paper. Section 3 presents the hybrid cold-plasma closure for Yang-Mills plasmas: the 
resulting Hamiltonian structure is explained in detail along with its momentum map property. The 
whole theory is directly applied to magnetized fluids with multipolar order and the resulting dynamics 
is then compared to Isayev's equations [38]. Section 4 couples the hybrid closure of the Vlasov- Yang- 
Mills equation with the dynamics of the Yang-Mills field in order to produce the hybrid formulation 
of chromohydrodynamics. Then, the same hybrid closure is applied to the case of ferromagnetic hy- 
drodynamics with quadrupolar interactions. Section 5 introduces the hybrid closure for systems with 
broken symmetry, such as nematic liquid crystals. In Section 6, the hybrid formulation of liquid crystal 
dynamics is presented, after deriving a suitable expression for the free energy which is inspired by the 
Frank energy. Then, the hybrid model of liquid crystals is shown to recover Eringen's formulation of 
micropolar liquid crystals [ ]. Applications to cubatic liquid crystals are also proposed. This paper 
focuses on Hamiltonian dynamics, thereby discarding the effects of dissipation. A geometric treatment 
of the dissipative dynamics in hybrid models is presented in [34]. 

2 Hybrid closure of the Yang-Mills Vlasov equation 

One easily verifies that the hybrid cold-plasma closure of the YM Vlasov equation 

/(x, p, a, t) = p(x, a, t) 5(p - P(x, t)) (14) 

yields another Hamiltonian fluid system, which allows for a microscopic statistical treatment of the YM 
charge dynamics. Indeed, upon introducing the fluid momentum 

m(x, t) = J p/(x, p, a, t) dpdcr = P(x, t) J p(x, a, t) da (15) 

and by using the relation 

p(x,a,t) = J /(x,p,<r,t) d 3 p, (16) 
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direct substitution of the chain rule formula 



5F _5F SF 
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into the YM Vlasov bracket (8) yields the Poisson bracket for the independent variables (p, m) 
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where the Lie bracket [V, V']^ := (V • V)V — (V • V)V is minus the Jacobi-Lie bracket on the Lie 
algebra X(M 3 ) of vector fields on R 3 [46, 35]. The above Poisson structure produces the equations 
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which constitute the most fundamental hybrid micro-macroscopic model. It is interesting to notice 
that the variables (p, m) are moments of the Vlasov distributions that belong to two different moment 
hierarchies, thereby producing one more reason for the word 'hybrid'. This model provides a new 
philosophy of approaching order parameter dynamics. Indeed, the macroscopic order parameter field is no 
more simply transported along the fluid, but rather the fluid continuously transports the whole statistical 
information about micromotion. This information is entirely encoded in the probability distribution 
function p. This type of hybrid closure has recently been introduced in the context of Smoluchowski 
dynamics of ferromagnetic particles [34]. 

Remark 1 (Lie-Poisson structure of the hybrid closure) We notice that the Poisson bracket (18) 
is Lie Poisson on the dual of the semidirect-product Lie algebra 

X(R 3 )(D^(M 3 ,C 00 (o*)) 

where J-"(IR 3 , C°°(o*)) denotes the space of scalar functions taking values in the Poisson algebra C°°(o*). 
Notice that this Lie algebra underlies the semidirect-product Lie group 

Diff(M 3 )(s)J r (M 3 ,C 00 (o*)) 

where the Lie algebra J^M 3 , C°°(o*)) also arises from an appropriate Lie group J^IR 3 ,^), where Q C 
Diff(T*0) is a Lie group whose properties have been recently investigated in [IS]. The Lie group Q 
can be regarded as a subgroup o/Diff can (T*0), i.e. canonical transformations on T*0. (More precisely, 
Diff can (T*0) should be replaced here by the group of strict contact transformations on T*OxR, see [18]). 
This group F{E?,Q) plays a similar role as the gauge group J-"(M 3 ,0) in chromohydrodynamics [21]. See 
also the gauge theory treatment of complex fluids in [12]. However, in the present case Q C Diff can (T*0) 
is infinite- dimensional, as opposed to O, which is finite- dimensional. 



2.1 Momentum map properties and singular solutions 

Since it is well known that moments of the Vlasov equation are momentum maps [22, 37], it becomes 
natural to ask whether the moment quantity p(x, a) is a momentum map. In particular, we notice that 
the moment p(f) is the zeroth-order moment of the moment hierarchy 

A n (x,a,t) = J P ®"/d 3 p, 
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whose momentum map property arises as the dual of the Lie algebra homomorphism {/3 n } l— ^2p® n — I f3 n , 
where _l denotes tensor contraction. This homomorphism associates to each symmetric contravariant 
tensor field (/? n (x, a)) ll '" %n d x i x . . . d x % n the Hamiltonian function /i(x, p, a) = p® n _l (3 n . Setting n = 
determines the Lie algebra inclusion i : J^M 3 , C°°(o*)) ^ C°°(T*M 3 x o*), such that i(/3o) = fio- Since 
the latter is evidently a Lie algebra homomorphism, then the dual map i*(f ) = f /d 3 p = Ao(x, a) is an 
equivariant momentum map. 

In order to understand the physical meaning of the hybrid closure, it is helpful to notice some relevant 
solutions of the kinetic equation for p in (19). First, one simply verifies that the solution 



p(x, a, t) = -D(x, t) 5 I a 



G(x, t) 



yields the CHD equation for the macroscopic YM charge density G G C° 
single-particle solution 

p(x,<r,t) = 6(x-Q(t)) S(a-T(t)) 



I s , o)*, while the Klimontovich 
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(x,a) = (x(t),T(t)) 



(x,a)=(x( i ),T(t)) 

thereby returning the fluid particle trajectory Q(t) together with the micromotion for its YM charge 
T(t). Another relevant class of singular solutions is provided by the expression 



p(x, a,t) = J w(a, s, t) <5(x - Q(s, t)) d k s . 



(21) 



Here Q : S M 3 is a time dependent map (embedding, Q G £1^(5*, M 3 )) from a lower dimensional 
submanifold S C ]R 3 into the physical space M 3 . On the other hand, the quantity w is interpreted as a 
map w : 5 — )• Den(o*) from 5 into the infinite-dimensional space Den(o*) of probability densities on the 
dual Lie algebra o*. Direct substitution of the above solution ansatz yields the following equations 



gQM) 

dt 

dw{a, s, t) 
dt 




Q(s,t) 

dw(a, s, t) d SH 
da ' da 5p 



--Q(s,t) 



(22) 



For example, if dim S = 1, the above solutions describe the evolution of a charged filament which supports 
the statistical distribution w(a, s, t) for the dynamics of the YM charge along the filament itself. In the 
simplest case when dimS" = 0, the label s is lost in (21) and this solution reduces to a point particle 
with position Q(t), supporting a distribution w(a, t) for the dynamics of its own YM charge. The latter 
example yields a system of equations composed of the ODE for the position coordinate Q and the PDE 
for the probability density w and this is a unique feature of the hybrid closure, whose singular solutions 
may be different from the usual Klimontovich particle solution (20). 

Remark 2 (Singular solutions are momentum maps) Notice that the singular solution (21) arises 
as a momentum map 

T*Emb(S,R 3 ) x T(S,C°°(o*)y ^ -F(M 3 , C°°(o*))* 
where the Poisson manifold P = T*Emb(S', M 3 ) x J-(S, C°°(o*))* is endowed with the Poisson bracket 



{F,K} 



fSF 5K 5K SF 



w { a 



d 5F d 5K 
da 5w ' da 5w 



ds da . 



(23) 
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The infinitesimal generator of the associated action reads as 



5F j)_ 



Hp[F] = - 
which arises from the group action 



w(s,a)£(Q(s),a)da + J ~ Ud* dm a)/da a,^\dsda, 



(QW, PW, «;<*>) = (QW, P(°) + (< W (°),(Ve X p(-^))exp(ta>oQ (0) ) , Ad 



exp( t i) 
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where Ad* is the group coadjoint operator on the Lie group Q underlying the Lie algebra C°°(o*) ; such that 
C°°(o*) = T e Q . Similarly, the notation (Vg _1 ) g makes use of the tangent lift of the right multiplication 
in Q. See [18] for more details on how the group Q is constructed. It is interesting to notice that 
the above construction coincides with the momentum map construction in [36] (in the absence of fluid 
motion) for geodesic flows on semidirect-product Lie groups, with the only difference that the group Q is 
infinite- dimensional in the present case. 

2.2 The mass density 

Although the hybrid equations (19) are physically consistent, in practical situations it is convenient to 
allow for the mass density 

D(x,t) = J p(x,a,t) da (24) 
to be a dynamical variable. This is possible upon writing the cold-plasma solution as 



f(x,p,a,t) = p(x,a,t)5[ p 



m(x, t) 



(25) 



so that the variables (D, p, m) are treated as three different variables. In this case, substitution of the 
chain rule formula 

5F 5F 5F SF 



Sf 5D + 5p +P 5m 



(26) 



in (8) yields the Poisson bracket 

{F,K}(m,p) = /m- 
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The resulting equations of motion read as 
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(27) 



(28) 



(29) 



which are accompanied by the constraint (24). The latter constraint belongs to a geometric class of con- 
straints arising as the zero- level set of a certain momentum map [46, 35]. Gauss' law in electromagnetism 
is a celebrated example of this construction [48]. 
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It is easy to derive the Kelvin circulation dynamics associated to equation (28). To see this, we begin 
by considering the general relation [46, 35] 



d 
dt 



( A( W= Ad g(f + ad gg- 



which holds for any Lie group element gG^ and any dual vector p £ q* in the dual of the Lie algebra 
= TffS ' . Then, by the Lie-Poisson nature of the bracket (27), one can write equation (28) in the form 



d_ 
dt 



(r]*m) = rf 
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~dt 
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where the star symbol denotes pullback and rj G Diff(IR 3 ) is a diffeomorphism such that 5H/5m = 77077 . 
Upon dividing by Dq = n*D, taking the circulation over a closed loop Cq yields 



1 d 

c D o dt 



77* (m • dx) 



c-o A) 



p V -r— da ■ dx 
dp 



Thus, changing variables produces the following circulation dynamics: 
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p V — - do- • dx , 
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where C = r\ o Co is a loop moving with the fluid velocity u = 77 o ry -1 = SH/5m. When the Hamiltonian 
H depends only on the first-order moment G = J apda, then the above relation recovers the circulation 
dynamics previously found in [24] for the Lagrangian dynamics of perfect complex fluids. 
Moreover, notice that the equations (28)- (29) leave the following quantity invariant: 



C(D, p) = j D A(D~ 1 p) do- dx . 



for any scalar function A. More particularly, the above quantity is a Casimir invariant for the Poisson 
bracket (27). This is easily seen by evaluating the bracket 



{F,C} 



I 



P\ 

A' (a 



D ( a 



d 5F OA' 
da bp ' da 
'dp d 5F 
da" 1 da bp 

D-iA'^A- 
da ' da 5p 

dDA d 6F 



da da 5p 



d 3 xda - J p £ 5F/5m A'd 3 xda - j D £ 5 F/s m J (A - pD~ l A') dad 3 x 
d 3 xda + j A'(£ SF/5m p- D £ 5F/5m (pD' 1 ) - pD~ 1 £ SF/Sm D^j d^xda 

d 3 xda + J A'(£ 5F/Sm p- £ SF/Sm (pDD- 1 ) )d 3 xdo- 
d 3 x do- + j A' (£ S F/5m P ~ £sF/5m p) d 3 x da 



=0, 



which holds for any functional F(m,D,p). The prime notation above A' denotes total derivative of A 
with respect to its argument. Thus, the hybrid closure possesses a large family of Casimir functionals 
that can be used for nonlinear stability studies [33]. 

Remark 3 (Alternative form of the hybrid closure) The previous section showed how the mass 
density D can be treated as an extra dynamical variable. However, this is possible only because of the 
constraint (24). This constraint prevents the variables D and p from being mutually independent, which 
instead would be preferable in a dynamical model. This problem can be easily avoided upon introducing 
the variable 

p(x,a, t) 



(p(x,a, t) 



D(x,i) 
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which is a scalar function on physical space 1R 3 taking values in the distributions on the dual Lie algebra 
o*. Compatibility with the constraint (24) requires that these distributions are normalized to unity at 
each point, that is 

J tp(x, a, t) da = 1 , Vx£l 3 , VtGl. 

This is easily interpreted by saying that each fluid particle carries a probability distribution on o* such that 
the total probability is the always the same for all particles. All the previous considerations concerning 
the equations (19) also apply to the above dynamics without essential modifications. In particular, the 
existence of singular solutions is preserved by this change of variables. 



2.3 Hybrid hydrodynamics of spin nematic phases 

While the magnetization dynamics has been addressed a great attention in various contexts, from fer- 
rofluids to spin glasses, the emergence of other magnetic phases with multipolar order has posed relevant 
problems towards dynamical modeling. For example, magnets with multipolar interactions have been 
recently studied in [56], where different multipolar ferromagnetic phases were presented together with 
their phase diagram. In particular, magnets with quadrupolar order are known under the name of spin 
nematic phases [51], since their order parameter is similar in form to that adopted for nematic liquid 
crystals. 

An interesting attempt for the introduction of quadrupole order in the dynamics was carried out in 
the late 90 's by Isayev and collaborators [38, 39, 9], who developed a Poisson bracket approach to include 
the traceless symmetric tensor 

Q = (a ®cr - ||tr| 2 I 

as an extra dynamical variable. Here, I is the identity matrix and the angle brackets (•) denote an 
averaging process over the space of spin, so that <x £ M 3 is the microscopic spin variable. (No confusion 
with the pairing bracket (-, •) should arise). Strictly speaking, one has <x G su(2) ~ So(3); however, in this 
paper we shall use the hat map <7y = —eijk°~k to identify so (3) ~ M 3 . Isayev's theory enlarges the spin 
symmetry to consider coadjoint motion on the Lie group SU(3) of special unitary transformations. As we 
shall see, this approach does not capture essential dynamical features, since it returns trivial dynamics 
for Q. Also, this approach is not feasible when higher order interactions must be considered and the use 
of higher-order tensors becomes necessary. For example, higher order multipolar interactions have been 
recently studied in [56], where different multipolar ferromagnetic phases were presented together with 
their phase diagram. 

From the preceding sections, it is clear that multipolar correlations are naturally encoded by hybrid 
models, which replace the magnetization order parameter by the distribution on the space of spin. For 
example, in the simplest case when both disclinations and external magnetic fields are absent, one can 
apply the Poisson bracket (27) to the following Hamiltonian function 

H{m,p) = \ j\u\ 2 d?*da + ^- j \M(p)\ 2 d 3 x + ^Jtt{Q 2 {p)) d 3 x , (30) 

which is the natural extension of the spin nematic Hamiltonian in [ ] to the case of an incompressible 
fluid with magnetization Ai = f pa da. Here, x an d p are left as unspecified physical constants. Here, 
the average defining the tensor Q is calculated with respect to the distribution p, so that 

Q(p) = [ p{a) (tr®tr- ||<x| 2 I ) da . (31) 



The Hamiltonian (30) can be regarded as the total energy for a (incompressible) spin nematic liquid 
without disclinations. Notice that inserting (30) in the first equation of (19) yields Euler's fluid equation 

— + (u-V)u = -Vp 
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for the divergence-less velocity u and this is due to the absence of both defects and external fields (which 
would produce non- vanishing force terms). On the other hand, upon using the Hamiltonian (30) in the 
second equation of (19), one obtains the following kinetic equation on the spin space 



dp 

at 



+ u-Vp = g^--a x (x^M + 2fiQa) 

OCT 



(32) 



where the gyromagnetic factor g has been introduced. Notice that standard tensor-index computations 
yield the pure transport relation dt + (u • V) A4 = 0, so that the initial condition Aio = reproduces 
spin nematic behavior [1, 43]. The same magnetization dynamics also arises from Isayev's theory [38], 
upon using the total energy (30). Moreover, upon taking the second-order moment of (32), one obtains 



OQa/3 

dt 



+ U • VQa/3 



9_ 

X 



M,Q 



2gn Q lv {{apa v (T 



ill ^a/i^y 



(33) 



where we have introduced the hat map notation Aiy :- 

Although the matrix commutator on the right hand side of (33) also emerges in Isayev's corresponding 
equation [38] 



80S 
dt 



+ u 



9_ 

X 



M , Q {1) + g^ Q$ (Sjpe^a + 5 av e lP ^)M 



nil) 



a/3 



(34) 



the latter replaces the other terms by a combination of the second and first moments, which leads to 
the trivial dynamics dfQ^ + u • \7Q^ = for spin nematics with vanishing magnetization. On the 
other hand, equation (33) requires computing the evolution of the moment hierarchy and special closure 
procedures need to be invoked. The study of these closure procedures and their application to realistic 
magnetic media will be the subject of future work. Another immediate advantage of the present hybrid 
setting is that it easily allows to incorporate the effects of higher order multipolar interactions (like those 
recently studied in [56]), whose corresponding terms would appear in the expression of the Hamiltonian. 
Notice that the family of Casimirs C = J A(p) dxdu allows a nonlinear stability analysis, which will be 
treated elsewhere. 



3 Coupling to the Yang-Mills field 

The preceding sections introduced the hybrid closure of the Yang-Mills Vlasov equation and investigated 
its geometric momentum map properties. This was done without considering the effects of the Yang- Mills 
fields. On the other hand, the YM Vlasov kinetic description first arose in the context of quark-gluon 
plasmas [ ], in which the field effects must be considered. This section shows how the hybrid closure 
combines with YM field dynamics. While this has mainly a pedagogical value, the implications of this 
construction in condensed matter systems are diverse. Indeed, the MHD limit of the whole Yang-Mills 
Vlasov system possesses the same Poisson bracket construction producing spin glass dynamics [30] and 
liquid crystal theories [2 ]. Thus, motivated by this suggestive analogy, we shall proceed by formulating 
the hybrid theory of Yang-Mills plasmas and by applying its MHD limit to ferromagnetic hydrodynamics 
with quadrupolar order. 



3.1 Hybrid description of Yang-Mills plasmas 

The YM Vlasov system is defined on the space Den(T*M 3 x o*) x T*^ 1 ^ 3 ^), where ^(M 3 ^) is the 
space of magnetic vector potentials. The Poisson bracket is given by the sum of a Vlasov term and 
(minus) a canonical one, involving the field variables (A,E): 

{F,K} = JfS^, S ^d 3 xd 3 pda-{F,K}(A,E) (35) 
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where one defines [21] 



and 



bF bK 
JJ'Jf 



{F,K}(A,E) 



bF bK 
Sf'Jf 



+ 



cr, 



d SF d bK 
dIrJf'daJf 



5F bK 
6A' ~bE 



5F bK 
6E' SA 



d 3 x 



The minus sign in (35) is consistent with the definition of the electric part E of the field as (minus) the 
momentum variable conjugate to the vector potential A [21]. The YM Vlasov Hamiltonian is given by 
(2). 

The new cold-plasma closure (25) generates a system on 36* (R 3 ) x Den(R 3 x o*) x T*0 1 (M 3 , o), where 
the space X*(R 3 ) of one- form densities on physical space contains the fluid momentum variable. The 
Hamiltonian (2) becomes 



H 



m 

D 



(a, A) 



d 3 x da + j D U(D) d 3 x - j p (a, A ) d 3 x da 

+ (\d A A\ 2 + |E| 2 ) d 3 x - J <E, d A A ) d 3 x (36) 



where the internal energy D U(D) has been inserted again to account for pressure effects. The new 
Poisson bracket in terms of the variables (D, p, m, A, E) reads 



{F,K} 



m 



5F bK 
()rn ' 5m 



d 3 x + I p\a, 

x 

bK 



d 5F d 5K 
da 5p ' da bp 



d 3 xda - {F,K}(A,E) 



P ( £SF/Sm - ^SK/Sm ~ ! d"xd<T 



6 Z) h»- 
bp) 



D £ 



■■SF/Sm 



5K 
bD 



£ 



5K/6m 



SF\ 

Id) 



d J x. 



(37) 



where one recalls the constraint D = j pda. This Poisson bracket is evidently the sum of (27) and 
— {F, G} (A, E). The equations of motions follow naturally by inserting the Hamiltonian (36) in the 
Poisson bracket (37). 

So far, the present construction yields the dynamics of the macroscopic fluid momentum m(x, t) = 
fpf(x,p,a,t)d 3 pda, which is not the variable commonly used in plasma physics. Indeed, a central 
role is usually played by the equation of motion for the macroscopic fluid velocity. At the particle level, 
the velocity coordinate is given by the minimal coupling formula v = p — (a, A), while we recall that 
in the macroscopic fluid setting of CHD one writes u = D~ 1 (m — (G, A)), see equation (11). Although 
the Hamiltonian (36) suggests the velocity definition u = D^ 1 (m — fp (a, A) da) , a careful analysis 
requires to apply the minimal coupling formula at the microscopic particle level, and to define the 
averaged macroscopic fluid velocity only at a second stage. Then, the hybrid equations can be written in 
terms of the macroscopic fluid velocity, which is expressed by using moments of the Vlasov distribution 
function. 

In analogy with the preceding discussion, one starts with the Poisson bracket (35) and the YM Vlasov 
Hamiltonian (2). Then, in order to introduce minimal coupling, one utilizes the well known relation 



P - W, A) 



(38) 



between momentum and velocity coordinates. Thus, one expresses the system in terms of the distribution 
f(x.,v,a,t). The Poisson bracket becomes 



«WC'f},-™ A)+ /'('-( 



d bF bK d bK bF 
dvJJ ' 5E ~ dvJJ ' bE 



(39) 
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where we define [50] 



(SF SK\ f SF SK\ 
1*7' *7V 1*7' *7li 



+ ( ex, B 



5 <5F 

v 

r7 ' ' ' dv Sf ' 3a sf 



d SK d SF 
dv Sf ' da Sf 



Once minimal coupling has been introduced, we are ready to derive the hybrid model for Yang-Mills 
plasmas by adopting the hybrid Vlasov closure 



f(x,v,a,t) = p(x,a,t) S[ v 



M(x,t) 
D(x,t) 



Upon inserting the chain rule relation 



5F 5F 5F 5F 

— = 1 h v ■ 

Sf 6D dp SM 



(40) 



(41) 



in (39), we obtain the hybrid Poisson bracket 



{ F,K } =f 



M 



SF SK 
SM' SM 



d s x 



D { £sf/6M - £&K/8M T~R ) ( "W 



5F_ 

Id 



j P (£sf/5M ^ - £sk/sm ^) d 3 xdcr - {F, K}(E, A) 



+ lp(tT, 



+ P{<T- 



d 5F d 5K 
da 8p ' da 5p 
5F d 5K 
SM'daJp 
5F 5K 5K 5F 
SM ' JE ~ SM ' SE 



A- 



d 3 xdcr + p(a,B 



A- 



5F 



SK d SF 
SM'daJp 

da dx 



SM SM/ 



da dx 



da 



(42) 



In the multi-fluid case, the above bracket represents the hybrid Yang-Mills version of the Spencer- 
Kaufman bracket for multicomponent magnetized plasmas [53, 54]. 
The Hamiltonian (2) now takes the hybrid form 



H 



1 



2 

I D~ J p{a ' A ° ) + J DU ^ + \ j (K A | 2 + I E | 2 )+ / (V,d A A c 



(43) 



so that Hamilton's equations become simply 

dD 

+ div(L>u) = 



du 



+ u • Vu 



dt 
D 

dp 
di 

which are to be accompanied by GauB Law 

SH 



dt 

+ div(pu) = (a, 



-Vp 

dp 
da' 7 



p a da, E + u x B 
A + A • u 



6 An 



= div A E+ pa da 
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and the field equations 

j p a da = curl A (ii A A)+ ad^ Q E - u J p a da 



<9E _ 5H bH 
~dt ~ IX ~ 5M 

bH a 

-8t = -jE = - E -" A ° 



The first moment of the p— equation yields the continuity equation in the covariant form 

d 



Of 



+ &d* Ao ) I pa da + div A ( u / p a da ) = . 



Notice that these equations contain precisely the same amount of physical information as the chromohy- 
drodynamics equations (11), (12) and (13) (see also [21]) because the dynamics of the Yang-Mills charge 
density G = J pa da does not involve any higher-order moment of p. 



Remark 4 (Hybrid Yang-Mills MHD) The MHD limit of chromohydrodynamics (Yang-Mills MHD, 
or YM-MHD) is known to possess many condensed matter applications. For example, the equations of 
spin glasses were shown [30] to possess the same Hamiltonian structure as Yang-Mills MHD. Thus, one 
is interested in finding the hybrid analogue of YM-MHD. This task can be easily accomplished by fol- 
lowing precisely the same steps as in [30], upon recalling that the total YM charge is now written as 
G(p) = Jpada. YM-MHD shares many analogies with ferromagnetic hydrodynamics [ ], which is 
generalized in the next section to systems with quadrupolar order. 



3.2 Ferromagnetic fluids with quadrupolar order 



The equations of ferromagnetic hydrodynamics (FHD) possess a well known Hamiltonian structure, 
which was discovered and investigated in [28]. More recently, a Vlasov equation of the same type as 
in [21] has been presented in relation to magnetized plasmas [5, 45]. In this section, we present an 
intermediate hybrid formulation of ferromagnetic continua, which illustrates a direct application of the 
hybrid approach to the case of quadrupolar interactions. For the purpose of simplicity, we shall consider 
incompressible fluid flows and the fluid momentum m will be denoted by u. 

Based on the work of Holm and Kupershmidt [28], it is easy to write the hybrid total energy for 
ferromagnetic (incompressible) hydrodynamics with quadrupolar order (cf. [1]): 



H 



lul 2 d 3 x 



1 

8^ 



J |B-4vrM(p)| 2 d ( jd 3 x+|yTr(Q 2 (p))d 3 x, 



(44) 



where the tensor Q is defined in (31) and a G so(3) ~ M 3 is the particle magnetic spin interacting with 
the magnetic induction field B. Evidently, when p = 0, the above expression recovers the incompressible 
version of the fluid Hamiltonian in [28]. Upon restricting to incompressible flows, the hybrid equations 
of motion arise from the Poisson bracket 



{F,K}(u,p,B) 



J5F SK\. 3 
u • curl 1 



d 5F d 5K 
da bp ' dcr bp 



d 3 xder 



f ( p SK P 5F \aZ a 

J P \^8F/Su - ^SK/Su I d Xdcr 

5K / 5 »bB/ 



B,£ 



bG 



<5F/<5u 



bB 



f 6F \d^ 



(45) 



which is derived from the Hamiltonian structure in [21] upon replacing = J a pd?a. For the Lie 
algebra so (3) ~ M 3 , the Lie bracket [•, •] above is given by the cross product, i.e. [£, u] = £ x v. Also, all 
the angle brackets (•, •) above are simply dot products. 



15 



Remark 5 (Ferromagnetic hydrodynamics v.s. Yang-Mills MHD) Notice that the pairings (■, ■} 
appearing in (45) allow for the more general case when the magnetic induction B is replaced by a Yang- 
Mills field taking values in the Lia algebra o, whose dual o* contains the Yang-Mills charge a. This 
produces precisely the Poisson bracket of hybrid YM-MHD. However, for ferromagnetic hydrodynamics, 
one has a G so(3)* ~ M 3 while the magnetic field is simply a two-form B-dS 6 $7 2 (M 3 ) (here, dS denotes 
an oriented surface element in so that the angle brackets (•, •) in (45) reduce to dot product. 



The equations of motion read 

— + (u-V)u + Vp = - / pV— da - B x curl — 
at J op oa 

dp _ dp d SH 

— + u-Vp = -ga- — x -a—j- 
ot oa oa dp 

<9B 

— + Vx(Bxu) = 

where the constant gyromagnetic ratio g has been inserted in the /5-equation for consistency. Upon using 
the hybrid Hamiltonian (44), redefining the pressure appropriately yields the hybrid equations 



+ (u • V)u + Vp = VB • J pa da - ^-B x V x(^B-Att j padc 



du 

M 

^ + div(pu) = g (B - AttM - 2pQa) -<tx^ 
ot ' oa 

<9B 

— + Vx(Bxu) = 
ot 

which evidently reduce to the incompressible form of the fluid equations in [28] by simply setting p = 0. 
We recall that the quantity C{p) = J A(p) dcrdx determines a whole family of Casimir invariants for the 
corresponding Poisson bracket (45). Again, these can be used to apply the energy-Casimir method by 
following [33]. Also, the magnetic helicity J* A • Bd 3 x (where B = curl A) is another Casimir for the 
Poisson structure of incompressible FHD, when the latter is expressed in terms of the vector potential 
A • dx G il 1 (M 3 ). While stability is left for future work, the next section extends the previous models to 
account for systems with symmetry breaking. More particularly, the case of liquid crystal flows will be 
considered. 



4 Hybrid models for systems with symmetry breaking 

In soft matter systems one often deals with broken symmetries. As briefly mentioned in the Introduction, 
this happens when the single particle dynamics is not entirely symmetric under its underlying configu- 
ration Lie group O. Usually the particle dynamics allows for a residual symmetry V, which is normally 
used to describe the particle micromotion on the symmetry-reduced phase space T*0/V ~ o* x (O/V). 
See e.g. [17, 24] for how the coset space emerges from the reduction of systems with symmetry break- 
ing. Here, the coset space M = O/V is the order parameter space, while O is sometimes called the 
broken symmetry group. The reduction theory for systems with broken symmetry has been developed in 
[15, 17, 24], to which we address for further details. In the remaining sections, the previous results on 
hybrid models will be extended to systems with broken symmetry. A hybrid formulation of liquid crystal 
dynamics will be presented and compared to Eringen's micropolar theory. Applications to cubatic liquid 
crystals with multipolar order will be proposed. 

4.1 Yang-Mills- Vlasov equation with symmetry breaking 

The YM Vlasov equation for systems with broken symmetry is written immediately by simply replacing 
the Lie- Poisson bracket on o* in (1) by the natural generalized Lie-Poisson bracket on o* x M [17]. This 
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leads to the following form of the YM Vlasov equation 



d_l 
dt 



+ /, 



SH_ 



+ cr 



df_ d_SH_ 
da' da Sf 



df fd_5H\ 
dn\daSfJ 



n + 



o_sh /an 

dn Sf \da J 



71 







(46) 



where n 6 M is the order parameter coordinate and the notation £m n stands for the infinitesimal action 
of O on M. Here the angle bracket (-, •) denotes either the pairing between the Lie algebra and its dual 
and the pairing between vectors and covectors on M. 

In order to obtain the fluid description of complex fluids with broken symmetry, one considers the 
following Vlasov distribution function /(x, p, n, a, t): 



/(x, p, a, n, t) = £>(x, t) S(n - n(x, £)) S( p 



m(x, t) 
D(x,t) 



Si a 



g(x,t) 



In a typical situation one has M = S 2 Cl 3 , that is the order parameter space M is the unit sphere S 2 
and points n £ S 2 on the sphere are identified with unit vectors in M 3 , so that £m n = £ x n. The fluid 
equations for an arbitrary Hamiltonian H(m,G, r\,D) read as [24, 15, 39] 



777 + £ 5H_ 

at 5m, 

777 ~r °£ 5H 
at 5m 

777 + I 5H_ ) rn 

at 5m, 



D = 



ad^H G + n o — 
jg on 



M 



■DV- 



51? 



G,V 



+ 



7n~ 



, Vn 



(47) 



where one defines the diamond operator o as (n o k, := (k, £m n )i f° r an Y covector k G T*M. Evidently, 
the quantity no k coincides with the momentum map for cotangent-lifted actions [46, 35]. The Poisson 
bracket corresponding to (47) is 



{F,K}(m,D,G,n)= J , 



I m • 


'SF 


SK~ 




f(c, 


'SF 


SK~ 








d 3 x + 








Sui' 


Sm 




SG 7 


Jg 



d 3 x 



D £ 



SK 

5F/5m ~ 

r f 5K 

&5F/5m 



£j\rf/Am 77T ) d 3 X 



£ 



SF_ 

SK/5m JJ^ 
SF 



5K/5m 



SG 



d 3 X 



+ 



SF 

5r7' 



SK\ 



SG J 



£sj< ) n 

5m 



/SK 
\ JrT' 



(-) 
\SGj_ 



M 



£sf_ ) n 

5m 



d 3 x, (48) 



This framework is particularly suitable for liquid crystals (the director field notation is n = n(x,t)), 
which will be the guiding example in this last section. Convenient Poisson bracket approaches to several 
types of liquid crystal dynamics were developed in [55, 12, 60, 39]. For example, the dynamics of uniaxial 
nematic liquid crystals is governed by the equations (47) and the Hamiltonian 



H= \ /^ d3x+ 2j / l G ! 2d3x + J DF(D-\ n,Vn)d 3 x. 



(49) 



where 



D F(D~ 1 ,n, Vn) = i^i(divn) 2 + (n • V x n) (iT 2 (n- V x n) + K' 2 ) + K 3 \n x V x n| 2 (50) 
is the Frank energy [7, 20]. 
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The presence in (50) of gradients of the order parameter field reflects the elastic properties of the 
liquid crystal material. In more generality, the emergence of these gradients is responsible for the failure 
of kinetic mean field approaches to complex fluid dynamics. Indeed, while mean field approaches are 
well known to apply to equilibrium situations, their employment in dynamical problems present several 
difficulties. An interesting attempt to overcome these difficulties is given by the following mesoscopic 
closure developed in [3, 4]: 

ft A ( A x( M(x,n,t) \ ( E(x,n,t) \ 
/(x, p, a, n, t) = p(x, n, t) 5 p — )5 [a — , 

where all hydrodynamical variables are defined on the extended configuration space M 3 x M. Thus, 
the microscopic information on the order parameter space M is retained by all dynamical variables. 
However, the resulting equations present several difficulties due to the high dimensionality of the space 
M 3 x M, which makes the dynamics difficult to solve even by using advanced numerical methods. 

The next sections show how these difficulties can be overcome by hybrid models, in which the 
macroscopic order parameter field is replaced by a position-dependent statistical distribution on the 
order parameter space M. 



4.2 Hybrid closure for systems with broken symmetry 

Recent computer simulations of cubatic liquid crystal phases [11, 2] and other complex fluids with 
multipolar order raise the necessity of accounting for the effects of correlations higher than second order. 
Previous sections showed how, for magnetic media, these correlation effects are easily taken into account 
by the hybrid approach. The following sections address the problem of extending this approach to 
systems with broken symmetry, with special emphasis on liquid crystal dynamics. 

A hybrid closure of the YM Vlasov equation for /(x, p,cr, n) can be easily obtained upon treating 
both coordinates (a, n) at the same level, thereby obtaining the closure 

, , . „/ m(x,t) 

/(x, p, a, n, t) = p(x, cr, n, t) 5 p 



D(x,t) 



with D(x) = J p(x, a, n) d<rdn. The above closure applies to nematic systems in which different molecules 
at the same point possess different orientation and angular momentum. This is a rather general situation. 
However, the position-dependent distribution /j(x, a, n, t) still depends on a high number of variables and 
thus it is not easy to handle. Thus, although the above hybrid closure may apply to nematic liquid crystal 
dynamics, it is more convenient to make simplifying assumptions. 

In the search for a hybrid description of complex fluids with symmetry breaking, one can seek an 
equation for a probability distribution on the space M 3 x M. This is precisely the same approach that is 
followed, for example, in ordinary Smoluchowsky theories, where such a distribution obeys a diffusion- 
type kinetic equation in the overdamped limit [10]. On the other hand, in this paper we consider 
systems with non-negligible inertial effects, so that the following results differ from the predictions of 
the Smoluchowsky approach. In the present setting, we shall consider the cold plasma closure 

/(x, P , a, n, t) = D(x, t) p(x, n,t)s(p- ^4) * ( ° ^ ' ' 



D(x, t) J V D(x,t) 

where J <p(n, n, t) dn = 1. This closure applies to situations in which different molecules at the same point 
possess different orientations and the same angular momentum. Thus, at each point in space one can 
have several rod-like molecules rotating in a moving plane with different orientations at all times. The 
resulting hybrid model possesses a Hamiltonian structure, which can be obtained by the same methods 
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that were presented in the first half of this paper. One obtains the following Poisson bracket 
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In the incompressible case (D = 1) the fluid momentum equals the velocity field u, i.e. m = u. Then, 
the hybrid equations of motion are 
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Notice that, when M = S 2 so that n is denoted by the unit vector n and the Yang-Mills charge G by G, 
the relation G • n = G • J rupdn = is not preserved by the hybrid dynamics, when the free energy J- 
contains moments of order higher than one. Thus, in this case, molecules that are initially uniaxial (i.e. 
satisfying G • Jn ip dn = 0) can change their configuration as time goes by. This is a typical situation 
in liquid crystal dynamics. For example, phase transitions (e.g. from uniaxial to biaxial) may occur in 
presence of defects, such as disclination lines. 

Also, notice that moments Jn . . .ntp dn of the distribution ip can be taken at any order, thereby 
generating a hierarchy that can be truncated at pleasure by simply assuming that the energy does not 
involve moments higher than a desired order. Indeed, the equation of the k-th moment does not involve 
higher order moments, as it usually happens in kinetic theory. Rather, the equation of the fc-th moment 
involves only moments of the same order k. Thus, the hybrid method proposed here always gives a 
closed system of equations for an arbitrary multipolar order parameter. Of particular interest in liquid 
crystal dynamics are moment quantities of the type 



Q(x, t) = J (I — a nn) tp(x, n, t) dn . 



(52) 



where J and a are constant factors. For example, when a = 1, the moment Q is known as microinertia 
tensor in Eringen's micropolar theory [13], while a = 3 yields a traceless alignment tensor that is used 
in the Landau-de Gennes theory [19]. It is perhaps not surprising that the equation for Q(x, t) turns 
out to be identical to that recently derived in [] 7] on a purely geometric basis. The same dynamics for 
Q was first found by Volovik in [60]. Indeed, the standard methods presented at the beginning of this 
paper can be used to show that the whole Poisson bracket (51) can be written in terms of the alignment 
tensor Q, thereby obtaining exactly the same Hamiltonian structure as in section 5.1.2 of [I 7]. Then, 
the expression of the free energy F(Q, VQ) usually coincides with the Landau-deGennes free energy [20]. 

The dynamics of Q is particularly important in the case when phase transitions (e.g. from uniaxial 
to biaxial) may occur in the system. This situation typically occurs in presence of disclination lines, 
which is an interesting case to study. This is the topic of the next section. 
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5 Hybrid description of liquid crystals with disclinations 



When orientational defects are present in liquid crystals, the dynamics of these defects is usually related 
to an equation of motion for a Yang-Mills magnetic field, which is known as disclination density. This is 
the well known gauge-theory approach to defect dynamics [40, 12, 59]. In the literature, the dynamics 
of the disclination density is given through the dynamics of the corresponding vector potential. In 
more geometric terms, the vector potential is a principal connection, while its covariant derivative is the 
curvature of the connection. Thus, the starting point is the emergence of a magnetic potential in the 
fluid model. This can be easily recognized by invoking reduction theory principles [24, 15] to write the 
field n G C°°(IR 3 , M) in terms of the order parameter group element x £ C°°(M 3 , O): 

where no is a (space-independent) fixed reference configuration, typically no = (0,0,1), and ??* is the 
push- forward of the diffeomorphism r/ = exp(5H/5m) transporting all fluid particles. Then, upon using 
standard properties of the push- forward [46, 35], simple differentiation yields 

Vn = ??* (Vx n ) = rj* {(^XX^ 1 ) %n ) = i]* (VxX^ 1 ) V*X"o =■ 7Af n , (53) 

where 7 = f]*{(^x) X 1 ) G r2 1 (M 3 ,o) is a magnetic vector potential. In the general case, one assumes 
that there exists a vector potential 7 such that the relation Vn = "y^n is verified at all times. As a 
result, the gradients of the order parameter field are naturally encoded in the vector potential, which is 
a new dynamical variable. 

5.1 Derivation of the Hamiltonian structure 

In the search of a hybrid model describing liquid crystal dynamics in presence of disclination lines, one 
can proceed by allowing for both magnetic and electric Yang-Mills fields and then operating on the 
Poisson bracket. One can start by considering the direct-sum bracket 

{F, K}(m, D, G, ^ A, E) = - {F, K] (A, E) + {F, K}(m, D, G, if) , 

where the second term on the right hand side is given by (51). Then, one can perform the change of 
variables 

G 1 — y G-div A E=:Q 

m i-)- m + E a x curl A a — A a divE a =: M . 

The momentum shift in the second line was already used for electromagnetic charged fluids by Holm 
[25, 26], who showed how this transformation leads to the Poisson bracket of plasma MHD. The first 
line is also used here to account for charge neutrality, i.e. div A E = 0. Both the above change of 
variables are given as shifts by a momentum map. Indeed, while the term — div A E is a momentum 
map arising from the action of the gauge group (as it appears also in electromagnetic plasmas [48]), 
the quantity E a x curl A a — A a div E a is a momentum map arising from the diffeomorphism action on 
the field variables. The same quantity also emerges naturally in electromagnetic charged fluids [25, 26]. 
This construction allows to apply a general result from [42] (see proposition 2.2), which generates the 
resulting Poisson bracket. Then, upon using the neutrality hypothesis 

div A E = => Q = G , 
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the assumption SF/S~Ei = 5K/5E = yields the final Poisson bracket for the MHD limit: 



{F,K}(m,D,G,p,j)= J 



m 



SF SK 
(5m' 5m 



dq + / (G, 



5F SK 



d 3 x 



nff 6K p 6F \a 

u I ^SF/Sm JJj - ^SK/Sm JJj I d< 7 



SK 

Jg 

SK 



£ 



SF 

8K/Sm 

5F\ 



d 3 x 



£sK/8mj^ J d 3 xdn 



SF 



<P \£ S F/Sm gtp 

[/ n SK 

J ill ^8F/8m - ^SK/Sm g 

d_SF fSK\ \ _ I d_5K_ 

dn S(p\ 5G J M I \dn Sp 



d 3 x 



SF\ 
6G )m 



d 3 xdn 



5F SK 

Jg'J^ 



(54) 



where we have replaced M by m for simplicity and A by 7 to comply with the notation used in the 
literature on liquid crystals. The resulting equations of motion read 
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(55) 
(56) 
(57) 
(58) 



Here the indexes a,b,c, . . . are Lie algebra indexes. These equations hold for an arbitrary form of the 
Hamiltonian H. Indeed, no specific assumption has been made on the form of the total energy. In the 
case of liquid crystals, one would hope to use Frank's expression of the free energy, which however still 
needs to be expressed in terms of the connection 7. The explicit form of Frank's energy is presented in 
the next section. 



Remark 6 (Disclinations) Disclinations are defined as discontinuities of the order parameter field n 
and the main difficulty involved in the study of their dynamics is that the Frank energy (50) diverges in the 
presence of discontinuities. The same problem emerges in the hybrid approach, because of the emergence 
of the gradient Vip, which blows up for (spatially) discontinuous distributions <p. This problem can be 
avoided by noticing that the relation 

is preserved by the dynamics, as it can be shown by a tedious computation ( or by general principles in 
reduction theory, see [l ] and remark 10). Then, if one replaces gradients according to the above 
relation, all divergences that appeared in the model are eliminated, as long as all other quantities are 
smooth. In [16], this fact was used to explore the various relations occurring among different dynamical 
models for liquid crystals. 

Remark 7 (Kelvin-Noether theorem) Notice that considerations analogue to those in Section 2.2 
lead the Kelvin circulation dynamics associated to equation (55). Indeed, upon repeating analogous steps, 
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one readily obtains the following Kelvin circulation dynamics: 

where j is a loop moving with the fluid velocity u = 5H/5m. 
5.2 The expression of the free energy 

In order to use the dynamical model in the previous section, it is necessary to express Frank's energy in 
the correct set of variables, including the connection describing defect dynamics. In Eringen's theory of 
liquid crystals [13], this quantity is denoted by 7 and it is given by jj = djXX T > due to the orthogonal 
property of x £ C°°(]R 3 , SO (3)), see equation (53). In order to express Frank's energy (50) in terms of 
7, one can use the relation 

djti = jj n = 7 ■ x n 

where the second equality makes convenient use of the well known hat-map isomorphism so (3) ~ R 3 
[46, 35], so that jj = 7 -. Upon fixing an appropriate basis {e a } on so (3) ~ M 3 , the connection form 
7 £ 17 1 (IR 3 ) <S> M 3 can be expressed as 7 = 7? e a dx l . In this context, the bold notation 7 that is ordinarily 
used for vectors in R 3 may become confusing and thus we shall relegate this notation to denote only 
componentwise quantities of the type 7^ = 7" e a or 7" = 7? dx l . For example, one computes 

n • V x n = n c e cjd djn d 

= n c£cjd£dablajnb 

= - n c n h (e dcj e dah )j a j 

= - n c n b (S ca 5 jb - 5 cb Sj a ) j a j 

= ~ naUjJaj + jl 

= — Tr(nn7) + Tr7 . 

By proceeding analogously, one has 

(n x V x n) i = e ijk nje kab d a n b 

= ^kij^kab n j^bcd1ca n d 

= - {0~ia0~jb - 0~ibO~ja) nj€ bcd J ca n d 

= ~ (ejcdjci + ticdjcj) n d rij 
= - e icd n d j cj nj 
= (nx {ijUj))., 

and also 

(divn) 2 = (diUi) 2 = (e iab j ai n b ) 2 = (7 • n) 2 , 

where we have defined the operator ffi := e^r/jfe, for any 3x3 matrix rjjk- Then, since 

|n x (jjUj) | 2 = \~fjTij f — (n • ("fjnj)) 2 = Tr(nn7 T 7) — Tr 2 (nn7) , 

we obtain the following alternative expression for the Frank free energy: 

F(n,7) = Ki (7 • n) 2 + (Tr(nn7) - Tr 7) (K 2 (TV(nn7) - Tr7) + K' 2 ) 

+ K 3 (Tr(nn7 T 7) - Tr 2 (nn7)) (59) 
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Therefore, a possible hybrid (quadratic) expression for the Frank free energy is given by 

T{ip, 7) = K\ J J <£>(x, n)(p(x, n') (7 • n) (7 ■ n') (n • n') dn' dn 

+ K 2 (p{x,n)(p(x,n')(Tr(nnj) - Tr 7) (Tr(n'n'7) -T^dn'dn 



+ 2-^2 / / ^( x > ™M X 5 n ') {n-n') (Tr(nn'7) + Tr(n'n7) - 2 (n • n') Tr7)dn'dn 



+ K 3 if(x, n)<p(x, ri) ( (n • n') Tr (nnV 7) - Tr(nn'7) Tr( n'nj) I dn' dn , (60) 



which evidently reduces to (59) upon restricting to the singular expression (p = 5(n — n(x)) (with n(x) £ 
S 2 ). The Mauer-Saupe factor n-n' has been introduced to make expression (60) symmetric with respect 
to permutations n o n' and Z2— reflections n — > —n (or n' — > —n'). Notice how the free energy (60) is 
completely determined by 7 and the quadrupole moment field J nmpdn. However, other forms of the 
hybrid free energy may also involve higher-order moments, similarly to what happens in the derivation 
of the celebrated Doi theory [10]. 

Remark 8 (Cubatic liquid crystals) Multipolar expressions emerge also in the description of cubatic 
liquid crystal phases [58, 11, 2], although expressing their free energy remains a difficult task. For 
example, in [11] the order parameter of cubatic phases is derived by making use of the ^-tensor 

35 f 1 
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whose dynamics can be derived from the hybrid approach by simply taking moments of the ^-equation, 
upon adding the appropriate terms to the expression of the free energy F(<p, 7). 

At this point, one would like to use the energy (60) to formulate an explicit model for the dynamics 
of liquid crystals in the presence of disclinations. However, the presence of defects is responsible for 
molecular configuration changes, thereby producing a varying microinertia tensor of the type (52), 
which modifies the expression of the total energy. 



5.3 Hybrid equations and relations to other models 

Upon restricting to the incompressible case, a possible hybrid expression of the Hamiltonian for liquid 
crystal dynamics is given by the functional 

H= X - y'|u(x)| 2 dx+^yG(x)-Q(^)- 1 G(x)dndx + Jt(<p, 7) dndx , (61) 

where the microinertia tensor Q(y) = J J <p (I — nn) dn evidently appears in the rotational energy, in 
analogy with rigid body dynamics. Notice that the above Hamiltonian assumes that the microinertia 
tensor is invertible at all times, while this is not necessarily true in the general case. At this point, 
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inserting the above Hamiltonian in equations (55)-(58) yields 

_ + (u . V )u-Vp = y — V^dn + 7 a xcurl— _ 7 "div — , 

^ + (u • V)G iGx (C^G) = (CT 1 *? • (I " nn)Q _1 G)n x || dn 

n x — — dn - div 7 — 
C/Ti dip 0*7 

g + »^ = (Q-. G ).„x* 

d t7 , + (di^) J:j + (u • V) 7 , = -diiQi^G) + 7j x (CP 1 *?) 

In order to get rid of the inverse matrix Q(y>) , we may introduce the angular velocity variable v in 
analogy with rigid body motion, such that G = Q((p)u. As a result, we have the final form of the hybrid 
equations: 



+ (u.V)u-Vp = -y— V^dn- 7 a xcurl— + 7 a div — , 



9u 

St 

^ + (u.V).)-QM^x I , = -^nxAgdn-div^ 



(62) 



|+u.V^ = ,.nx^, (64) 
d t Ji + (diU 3 )^ + (u • V)7 i = -div + 7i x i/ (65) 
which follow upon making use of the relations 



J j + (u • V)^ (I — nn) udn = v X Qi/ 



obtained from equation (64) by standard tensor index computations. 

At this point, one may wonder how equations (62)-(65) are related to ordinary models for complex 
fluids, under the substitution ip = 5(n — n(x, t)), with n(x, t) G S 2 . Indeed, since the Free energy was 
derived from the Frank energy upon using dixi = 7^ x n, the above hybrid model is compatible with 
Ericksen-Leslie dynamics [ ]. In particular, if one assumes that no defects are present at the initial 
time (i.e. 7(x, 0) = 0) and upon restricting to rod- like molecules (i.e. v ■ n = 0), this compatibility is 
ensured by geometric symmetry arguments [15, 16, 17, 24]. 

An interesting observation also arises when the hybrid equations are written in terms of the moments 
Jnn...n<p(x,n,t)dn. More particularly, since the free energy (60) depends only on second-order 
moments, then the total energy can be completely expressed in terms of the symmetric microinertia 
tensor (52) (with a = 1). In this way, the Poisson bracket (54) can be transformed into the Poisson 
bracket of Eringen's equations for micropolar media [13] (see section 8.11.2 in [ ]). Therefore, the 
hybrid equations of motion recover Eringen's description of micropolar liquid crystals. Notice that the 
same bracket structure involving the tensor (52) appeared previously in [60]. 

Higher moment orders may also appear in the expression of the free energy. Then, the moment 
hierarchy can be truncated at the highest order and the resulting theory generalizes ordinary fluid 
models of liquid crystal dynamics. Indeed, upon taking the k-th order moment of equation (64), the 
resulting moment equation does not involve any moment of order higher than k, differently from other 
usual kinetic theory approaches. This result may be used, for example, for defect dynamics in cubatic 
liquid crystals, whose order parameter can be defined as a combination of higher order moments, as 
shown in [11]. This interesting possibility will be explored elsewhere. 
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Remark 9 (Differences with Doi's theory) Besides considering dominant inertial effects and ne- 
glecting dissipation, the present approach differs from Doi's theory of polymeric fluids [10] also because it 
encodes disclination dynamics. This has always been an outstanding task in Smoluchow ski-like models, 
since it is not clear how disclinations can be incorporated into the framework of kinetic theory. The 
present hybrid approach solves this problem in the energy- conserving case, in which moment equations 
close naturally at any order. 

Remark 10 (The action of the order parameter gauge group) Notice that the construction above 
can be easily generalized to an arbitrary order parameter group O and it is geometrically equivalent to the 
usual geometric construction of complex fluid dynamics [24, 15], including the affine action of the gauge 
group J"(R 3 , O) on the space J7 1 (IR 3 ,o) of connection one-forms. Indeed, replacing the order parameter 
field n G ^(M 3 , M) by the distribution-valued function ip G J-(M. 3 , Den(M)) has the only result of changing 
the action of the gauge group J-(R. s ,0) on the respective quantity. For example, while J-(M?, SO(3)) acts 
on n(x) by matrix multiplication (i.e. n \-¥ %n ; with x £ -^(M 3 , SO(3))), the gauge action on ip(x,n) is 
<p(x,n) H> ipi^x^n). 

6 Conclusions 

This paper has introduced a new approach to complex fluids with multipolar order in which the fluid 
transports the statistical information that is encoded in a position-dependent probability distribution. 
In the context of Yang-Mills plasmas, this approach arises as a cold-plasma- like closure of the underlying 
Vlasov kinetic equation, which possesses momentum map properties. This closure involves moments 
belonging to different hierarchies thereby producing a hybrid description of the YM plasma. The resulting 
equations of motion possess a natural Hamiltonian structure, which in turn possesses a whole family of 
Casimir invariants. As a first application, the example of ferromagnetic hydrodynamics with quadrupolar 
interactions has been presented and compared with previous results in the literature [38]. 

The second part of the paper extended the hybrid formulation of complex fluids to systems exhibiting 
symmetry breaking. Again, starting from a microscopic kinetic approach, a hybrid moment closure has 
been presented, which involves a statistical distribution on the order parameter space. The case of 
broken symmetry presents some difficulties that where overcome by encoding the gradients of the order 
parameter field in a new variable that accounts for defect dynamics. Motivated by cubatic liquid crystal 
phases, the dynamics of liquid crystals was taken as the guiding example and a hybrid expression of the 
Frank free energy was presented. The resulting hybrid liquid crystal model is compatible with Ericksen- 
Leslie dynamics and it recovers Eringen's formulation of micropolar liquid crystals [13] when the total 
energy can be expressed in terms of second order moments. The same hybrid model also recovers 
equations 25 in [60]. The possibility of higher order moments appearing in the free energy generalizes 
the theory to account for the higher correlation effects that appear in cubatic phases [58]. 
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A A generalization of defects 

In many situations, the elastic energy of complex fluids depends on the gradients of elements O G 
J"(IR 3 ,e>) in the gauge group. After symmetry is taken into account, this leads to the introduction of a 
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connection one form 7(x) = _1 VO 6 Q 1 (IR 3 , o), which is also known as wryness tensor in the context 
of liquid crystal dynamics [13]. For example, when O = 50(3), the quantity 7 represents orientational 
defects, such as disclinations whose density is given by the associated curvature B = cP'j [12, 59]. A 
famous example of how this quantity appears in complex fluids is found in the dynamics of spin glasses 
[12, 59, 30]. 

Given the similarity with the gauge approach to complex fluids (see remark 1), the hybrid construction 
can be extended to include a connection taking values in the Poisson algebra C°°(o*). In this context, 
the gauge group can be considered to take values in the infinite dimensional (sub)group of canonical 
transformations Diff can (T*0), so that at every point x in physical space there is an associated canonical 
transformation g(x.) that is responsible for the statistical micromotion. Then, if gradients of the type 
Vg appear in the total energy of the system, one is led to the introduction of the one form 

0(x,a) = g^Vg e n x (R 3 ,C°°(o*)). 

This quantity is completely absent in soft matter modeling and emerges in the present context for the 
first time. Its dynamics can be easily derived by applying the affine reduction presented in [I 5], so that 
the Lie-Poisson bracket 
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where summation over repeated indexes has been used. In the absence of fluid flow (i.e. m = 0) and in 
one spatial dimension, the last two equations above belong to a class of systems whose integrable cases 
were recently studied in [27]. 

In the context of complex fluids, the connection one form ■& encodes defects associated to the phase 
space dynamics of the order parameter and its physical implications are yet unknown. Indeed, to the 
author's knowledge, no Hamiltonian H(m, p, t3) is available in the physics literature. As a conclusive 
remark, we notice that if such a Hamiltonian depends on p only through the first moment G = J a pda 
and ■& is linear in a (i.e. $(x, a) = (a, 7(x))), then the above dynamics reduces precisely to the 
Hamiltonian structure of Yang-Mills MHD [30, 31], with YM magnetic field B = d 7 7 (also known 
as disclination density). 
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